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We report the first six-dimensional quantum dynamical cal¬ 
culations of dissociative adsorption and associative desorp¬ 
tion. Using a potential energy surface obtained by density 
functional theory calculations, we show that the initial de¬ 
crease of the sticking probability with increasing kinetic en¬ 
ergy in the system H2/Pd(100), which is usually attributed 
to the existence of a molecular adsorption state, is due to 
dynamical steering. In addition, we examine the influence of 
rotational motion and orientation of the hydrogen molecule 
on adsorption and desorption. 


The dissociative adsorption and associative desorption 
of hydrogen on transition metals is of direct technological 
relevance as well as of general interest for the very fun¬ 
damental understanding of elemental chemical processes 
occurring in catalysis and surface chemistry. Molecu¬ 
lar beam experiments of the dissociative adsorption of 
H 2 on various transition metal surfaces like PdflOO) [lj, 
Pd(lll) and Pd(110) §, W(lll) [|, W(100) |,|,|, and 
Pt(100) Q revealed that the sticking probability initially 
decreases with increasing kinetic energy of the beam. 
Such a behavior is usually explained by a precursor mech¬ 
anism j7|]. In this concept the molecules are temporar¬ 
ily trapped in a molecular physisorption state, the so- 
called precursor state, before they dissociatively adsorb. 
The trapping occurs because the molecules lose energy 
when impinging on the surface, mainly to the substrate 
phonons |ij . The energy dependence of the sticking prob¬ 
ability is thus related to the trapping probability into the 
precursor state, and it is this trapping probability which 
decreases with increasing energy. In all experiments men¬ 
tioned above the sticking probability rises again at higher 
energies which is interpreted as such that adsorption via 
direct activated paths then becomes dominant. 

The relevance of the precursor mechanism for the disso¬ 
ciative adsorption has been discussed for a long time, and 
the discussion remains quite lively [|]j§]- One may argue 
that especially for hydrogen adsorption on metal surfaces 
the energy transfer to the phonon system should be rel¬ 
atively small due to the large mismatch in the masses 
of adsorbate and substrate making the precursor mecha¬ 
nism ineffective. Therefore also direct non-activated ad¬ 
sorption together with a dynamical steering effect has 
been suggested in order to explain the initial decrease of 
the sticking probability (see, e.g., ||||), however, there 


had been no theoretical analysis of the importance of 
this effect. Recent density functional theory calculations 
|JTo| of the potential energy surface (PES) of H 2 /Pd(100) 
show that there exist non-activated as well as activated 
paths to dissociative adsorption in this system, but they 
do not yield any molecular adsorption well. These re¬ 
sults motivated dynamical calculations using a model 
PES with activated as well as non-activated path to dis¬ 
sociative adsorption Q, in which three degrees of free¬ 
dom of the hydrogen molecule were taken into account. 
This investigation suggested that the steering mechanism 
could be operative in the system H2/Pd(100), but the 
results showed large quantitative discrepancies to the ex¬ 
periment [jiijl - 

In order to clarify the nature of the dissociation pro¬ 
cess, we parametrized the ab initio PES of H2/Pd(100) 
thereby obtaining a well-defined function describing the 
potential energy surface on which the hydrogen molecule 
moves. This function depends on six coordinates, and 
together with the kinetic energy operator it defines the 
Hamiltonian for the nuclear motion. The Hamiltonian 
is then used for a quantum dynamical study where for 
the first time all six degrees of the hydrogen molecule 
have been treated fully quantum mechanically in a simu¬ 
lation of dissociative adsorption and associative desorp¬ 
tion. The specific coordinates considered in our study are 
the three coordinates of the center of mass of the molecule 
X, Y , and Z (Z is the distance from the surface), the az¬ 
imuthal orientation </> and the polar orientation 9, and 
the H-H distance d. The potential in the Zd plane is de¬ 
scribed in reaction path coordinates s along the reaction 
path and r perpendicular to it ]l2| ] . The potential is then 
given in the following form 

V (X,Y,s,r,&,</>) = V corr + V rot + V vlb (1) 
with 
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G = 2tt /a is the basis vector of the reciprocal lattice, a 
is the lattice constant, /i is the reduced mass of the hy¬ 
drogen molecule, and oj(s) is the frequency perpendicular 
to the reaction path. The functions Vm?n(s), oj(s), and 
A r(X, Y, s) are determined such that the difference be¬ 
tween V(X,Y,s,r,8,<j)) and the ab initio total energies, 
which have been calculated for more than 250 configura¬ 
tions, on the average is smaller than 25 meV. 

Although the PES does not have a precursor molecular 
adsorption state, the results of the six-dimensional quan¬ 
tum dynamical calculations concerning the dependence 
of the sticking probability on the kinetic energy in the 
system H 2 /Pd(100) agree very well with the measured 
results. This will be explained by a purely dynamically 
steering effect. In addition we will show that rotations 
in general suppress sticking, but that steric effects, i. e., 
effects due to the orientation of the molecules, can lead 
to a strong enhancement of the sticking probability. 

The quantum dynamics is determined by solving the 
time-independent Schrodinger equation for the two hy¬ 
drogen nuclei moving on the six-dimensional PES. We 
use the concept of the local reflection matrix (LORE) 
[ 131 and the inverse local transmission matrix (INTRA) 
[ 141. This very stable coupled-channel method has been 
employed before in a high-dimensional study of the ad¬ 
sorption of H 2 /Cu( 111 ), where, however, the polar ori¬ 
entation of the molecule has not been varied [^5|. In 
the present calculation up to 21,000 channels per total 
energy are taken into account; with channels the eigen¬ 
functions of the asymptotic Hamiltonian are meant. The 
PES is obtained using density functional theory within 
the generalized gradient approximation (GGA) |^6| and 
the full-potential linear augmented plane wave method 
(FP-LAPW) (see Ref. juj and references therein). 

Figure ^ presents our results for the sticking probabil¬ 
ity as a function of the kinetic energy of the incident H 2 
beam. The dashed curve corresponds to H 2 molecules 
initially in the rotational ground state j. c = 0. This curve 
exhibits a strong oscillatory structure for low energies 
in contrast to the experimental results. Such oscillatory 
structures have already been found in three-dimensional 
model calculations using potential energy surfaces with 
activated as well as non-activated paths to adsorption 
|n],[r^]. For a PES without a molecular adsorption well 
the oscillations are related to quantum effects of the hy¬ 
drogen beam, namely the opening up of new scattering 
channels for increasing incident energy |11|. 

In molecular beam experiments the incident molecules 
are not all in the rotational ground state, the rotational 
states are occupied according to a Boltzmann distribu¬ 
tion with a temperature of 0.8 of the nozzle temper¬ 
ature Q. In addition, the beam is not strictly mo- 
noenergetic, but has a certain energy spread. Both ef¬ 
fects, the rotational distribution and the energy spread, 


are incorporated in the solid curve where we have as¬ 
sumed a typical energy spread of molecular beams of 
A E/Ei = 2A v/vi = 0.2 JlJ ( E, and Vi are the initial 
kinetic energy and velocity, respectively). This curve 
should be compared with the experimental results of Ren- 
dulic et al. 0. In view of the fact that no parameters are 
used in our calculations, the agreement between theory 
and experiment is quite satisfactory. The initial decrease 
of the sticking probability is well reproduced although 
no precursor state exists in our PES. Hence a purely dy¬ 
namical mechanism has to be the origin of this energy 
dependence. 

The mechanism, a dynamical steering effect, is dis¬ 
closed by the fact that more channels in the dynami¬ 
cal simulation are needed at low energies than at high 
energies. Usually it is the other way around in coupled- 
channel calculations since the higher the energy, the more 
channels become energetically accessible. Thus the cal¬ 
culations reflect that there is a strong rearrangement be¬ 
tween the different channels at low energies due to the 
steering. At high energies the dynamics is closer to the 
adiabatic limit which means that during the dissociation 
event there is little redistribution among the channels. 

An illustration of the steering effect is given in fig. 0 
There a two-dimensional cut through the six-dimensional 
PES is shown together with two typical trajectories of 
impinging H 2 molecules. Besides a difference in velocity, 
both trajectories have equivalent impact parameters. In 
the case of the left trajectory the incident molecule is so 
slow that the forces acting on the molecule can steer it 
towards a configuration where non-activated dissociative 
adsorption is possible, while in the case of the right tra¬ 
jectory the molecule is too fast for the forces to divert 
the molecule significantly, it hits the repulsive part of 
the potential and is reflected back into the gas phase. If 
the incident energy of the molecule is further increased, 
from a certain energy on it will have enough energy to 
directly traverse the barrier. This leads to the increase 
of the sticking probability at higher energies (see fig. [f). 

Of course in a complete description of the adsorption 
process also energy losses of the impinging molecules to 
the substrate phonons have to be considered, in spite of 
the large mass mismatch which makes the energy trans¬ 
fer rather small. Energy losses of the molecules would 
slow down the molecule and thus make the steering even 
more effective. Furthermore, the inclusion of the sub¬ 
strate motion would broaden the velocity spread of the 
incident molecules relative to the substrate atoms |T0. 

Figure [i] also demonstrates that by taking into ac¬ 
count the rotational population of the incoming beam in 
adsorption, the averaged sticking probability is slightly 
decreased as compared to molecules in the rotational 
ground state. This effect is shown in more detail in fig. 0 
which displays the sticking probability versus initial ro¬ 
tational quantum state for a fixed incident kinetic energy 
of Ei = 0.19 eV. The diamonds correspond to the orien- 
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which decreases with increasing ji. This means that the 
faster the molecules are rotating, the more the dissocia¬ 
tive adsorption is suppressed, and is caused by orienta¬ 
tional hindering : molecules with a high angular mo¬ 
mentum will rotate out of a favorable orientation towards 
adsorption during the dissociation event. 

From fig. and the principle of microscopic reversibil¬ 
ity, it follows that the population of rotational states in 
desorption should be lower than expected for molecules in 
thermal equilibrium with the surface temperature. This 
so-called rotational cooling has indeed been found for H 2 
molecules desorbing from Pd(100) Q. Figure || shows a 
comparison between theory and experiment of the rota¬ 
tional temperatures T ro t of desorbing molecules. These 
temperatures are determined by fitting the rotational 
state population in a Boltzmann plot to a straight line. 
The slope then gives (—fcsTrot)- In order to compare 
the calculations with the experimental data which were 
limited to low j states, for the theoretical analysis only 
the lowest four rotational states have been considered. 
There is again a rather good agreement between theory 
and experiment. 

The results of figs, [l] and [| have been determined by 
summing over all azimuthal quantum numbers to. From 
the PES we know, however, that the most favorable con¬ 
figuration towards dissociative adsorption is met when 
the molecular axis is parallel to the surface. Rotational 
eigenstates with azimuthal quantum number m = j have 
their rotational axis preferentially perpendicular to the 
surface, so that the molecular axis is oriented mainly 
parallel to the surface. Molecules rotating in this so- 
called helicopter fashion should dissociate more easily 
than molecules rotating in the cartwheel fashion with 
the rotational axis preferentially parallel to the surface 
(m = 0) since the latter have a high probability hitting 
the surface in an upright orientation in which they can¬ 
not dissociate. This steric effect has been investigated in 
a number of model studies for activated adsorption em¬ 
ploying three and four-dimensional quantum dynamics 
23}| or mixed classical-quantum dynamics |f2d|| . 
also the sticking probability for initial rota¬ 
tional quantum state for nii = 0 and m; = ji has been 
plotted. The results show that the steric effect is oper¬ 
ative in systems with activated as well as non-activated 
paths towards dissociative adsorption. Indeed the differ¬ 
ence in the sticking probability between initial mi = 0 
and to,; = ji rotations increases with increasing quan¬ 
tum number j t . This reflects the fact that the larger the 
rotational quantum number j is, the larger is the differ¬ 
ence in the orientation between to = 0 and to = j states. 
The to, = ji data demonstrate that the steric effect even 



over-compensates the suppression of the sticking prob¬ 
ability by rotational motion since the nii = ji sticking 
probability increases with increasing ji. 

In conclusion, we reported the first six-dimensional 
quantum dynamical study of dissociative adsorption 
on and associative desorption from surfaces. We em¬ 
ployed a potential energy surface obtained by detailed 
density functional theory calculations for the system 
H2/Pd(100). The results show an initial decrease of the 
sticking probability with increasing energy which, in con¬ 
trast to common believe, is here not due to a molecular 
adsorption state, but can be explained by a dynamical 
steering effect. Rotational motion in general suppresses 
sticking, however, a strong enhancement in the sticking 
probability for molecules with their axis parallel to the 
surface is found. The results clearly demonstrate that a 
high-dimensional dynamical analysis of the PES is indeed 
important and provides, due to the microscopic informa¬ 
tion, quantitative as well as new qualitative understanding 
of processes at surfaces. 
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FIG. 1. Sticking probability versus kinetic energy for a H 2 
beam under normal incidence on a Pd(100) surface. Dashed 
line: H 2 molecules initially in the rotational ground state; 
Solid line: H 2 molecules with an initial rotational and energy 
distribution adequate for molecular beam experiments (see 
text); circles: experiment (from ref. [jlj]). 


FIG. 2. Illustration of the steering effect. The con¬ 
tour lines show the potential energy surface in eV along a 
two-dimensional cut through the six-dimensional configura¬ 
tion space of a hydrogen molecule. The molecular axis is 
parallel to the surface. The Pd atoms are at the potential 
maxima. The surface coordinate corresponds to the X-axis, 
while the reaction path coordinate s for s —> 00 is equiva¬ 
lent to the H 2 center of mass distance from the surface, for 
s < —2.5 A to the distance between two adsorbed H atoms. 
The left trajectory corresponds to a slow molecule (kinetic en¬ 
ergy Ekin = 0.05 eV), the right trajectory to a fast molecule 
(Ekin = 0.15 eV). 


FIG. 3. Sticking probability versus initial rotational quan¬ 
tum state ji. Diamonds: orientationally averaged sticking 
probability (eq. |H|), triangles: rrii = 0 (cartwheel rotation), 
circles: rrii = ji (helicopter rotation). The initial kinetic en¬ 
ergy is Ei = 0.19 eV. 


FIG. 4. Rotational temperatures of desorbing H 2 
molecules. Circles: experiment (from ref. |20|), diamonds: 
theory. The solid line corresponds to molecules in equilib¬ 
rium with the surface temperature. 
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